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ABSTRACT 

There are in the hterature several theories to explain the mass loss in stellar 
winds. In particular, for late-type stars, some authors have proposed a wind 
model driven by an outward-directed flux of damped Alfven waves. The winds of 
these stars present great amounts of dust particles that, if charged, can give rise 
to new wave modes or modify the pre-existing ones. In this work, we study how 
the dust can affect the propagation of Alfven waves in these winds taking into 
account a speciflc damping mechanism, dust-cyclotron damping. This damping 
affects the Alfven wave propagation near the dust-cyclotron frequency. Hence, if 
we assume a dust size distribution, the damping occurs over a broad band of wave 
frequencies. In this work, we present a model of Alfven wave-driven winds using 
the dust-cyclotron damping mechanism. On the basis of coronal holes in the Sun, 
which present a superradial expansion, our model also assumes a diverging ge- 
ometry for the magnetic fleld. Thus, the mass, momentum, and energy equations 
are obtained and then solved in a self-consistent approach. Our results of wind 
velocity and temperature proflles for a typical K5 supergiant star shows compat- 
ibility with observations. We also show that, considering the presence of charged 
dust particles, the wave flux is less damped due to the dust-cyclotron damping 
than it would be if we consider some other damping mechanisms studied in the 
literature, such as nonlinear damping, resonant surface damping, and turbulent 
damping. 

Subject headings: MHD — plasmas — stars: winds, outflows — stars: late-type 
— supergiants — dust, extinction 
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1. INTRODUCTION 

Cool supergiant stars present high mass loss rates, M ~ 10^*^ — 10^^ M0 yr^-*^, and low 
terminal velocities, u^o ~ 10 — 100 km s^^ (Cassinelli 1979), which are less than the escape 
velocity at the stellar surface (feo)- Typically, it is observed that Mqo ^ Although 
the physical mechanism that drives these winds is still poorly known, driving the wind 
with magnetohydrodynamic waves is one of the most promising mechanism proposed to 
date. Several authors (e.g., Liberatore et al. 2001; Elitzur & Ivezic 2001; Hofner et al. 2003; 
Woitke & Niccolini 2005), motivated by the observed high stellar luminosity and low effective 
temperature of these stars, have proposed radiatively dust driven wind models to explain the 
wind acceleration. Acoustic waves generated at a pulsating phase could provide a propitious 
ambient medium for dust formation; this is, a high-density region with low temperature. 
Hence, radiative pressure acting on dust particles could accelerate the grain and, if the gas 
and the dust were dynamically coupled, grains could drag the gas and provide the mass loss. 
However, this mechanism fails when gas and dust are not well coupled (Sandin & Hofner 
2003) . Another faihng aspect of this mechanism is that there are observations attesting that 
the wind forms before the grain formation point (Carpenter et al. 1995). 

In this aspect, Alfven waves might play an important role in these scenarios, producing 
significant mass loss and low terminal velocities (Lamers & CassineUi 1999, p. 294). In this 
mechanism, the wave magnetic pressure gradient is added into the momentum equation in 
order to drive the wind. 

Hartmann & MacGregor (1980) investigated the effects of Alfven waves on the outer at- 
mospheres of cool supergiant stars. They estimated an initial wave flux of ~ 10^ ergs cm"^ s~^ 
This flux is of the same order of magnitude as the one estimated for the solar wind (Baner- 
jee et al. 1998). However, they found that undamped waves led to very high terminal wind 
velocities. In order to obtain terminal velocities that were in reasonable agreement with 
observations of cool supergiants stars, it was necessary to dissipate the waves on a length 
scale of ~ Tq {tq is the star radius). The physical mechanism responsible for the damping 
was left unspecified. Holzer et al. (1983) argued that this length scale was very restrictive 
and that we could not observe a large variety of stars if we had such selective dissipation 
scale for the waves. This problem was known as the fine tuning problem. Jatenco-Pereira & 
Opher (1989) showed that flux tubes expanding superradially from the stellar surface, as is 
observed to happen in solar coronal holes (Cranmer Sz van Ballegooijen 2005), could account 
for a broad band of dissipation scales, overcoming this problem. 

In Falceta-Gongalves & Jatenco-Pereira (2002, hereafter Paper I), a hybrid model was 
presented, in which the effects of a flux of Alfven waves acting together with radiation 
pressure on grains were analyzed as acceleration mechanisms of winds of cool supergiant stars. 



-3- 



In that model, an empirical temperature profile was introduced. The presence of grains had 
two roles in that model: grains were used to simulate a strong damping mechanism for the 
waves and as an additional acceleration mechanism due to the exertion of radiation pressure 
on them. Hence, a more realistic wind model showing results consistent with observations 
was obtained. 

Here, wc present a model in which a flux of Alfven waves propagating outward from the 
star is damped due to the interaction between the waves and charged dust particles, thus, 
accelerating the wind. In order to obtain both the temperature and the velocity profiles, 
we solve the mass, the momentum, and the energy equations of the wind. In this work, we 
improve the model presented in Paper 1 by: (1) inserting the energy equation into the model 
in order to obtain the temperature profile of the wind; and (2) inserting a more realistic 
damping mechanism for the waves due to the presence of grains. Here, we examine the 
infiuence of this damping mechanism instead of simulating a strong damping. In §2, we 
describe the model used to accelerate the wind: we describe the modified wave dispersion 
relation due to the presence of charged dust in §2.1, the geometry considered for the magnetic 
field lines in §2.2, the equations that describe the wind dynamics in §2.3, and the radiation 
pressure on dust in §2.4. In §3 and §4, we present our results and conclusions. 



2. MODEL 

2.1. Dust- Cyclotron Damping 

The propagation and damping of Alfven waves in a dusty plasma has been considered by 
many authors (Pillip et al. 1987; Mcndis & Rosenberg 1992; Shukla 1992; Falceta-Gongalves 
et al. 2003). Dust grains immersed in an ambient plasma become charged due to plasma ion 
and electron fluxes into grain surfaces. Although the number of dust particles is smaller than 
the number of ions, the process of dust charging is efficient and these particles can obtain 
charges qd = —ZdC for Zd on the order of ~ 10° — 10^ in astrophysical media (Goertz 1989; 
Mendis & Rosenberg 1994). Once charged, these particles suffer the infiuence of the magnetic 
field which gives rise to a cyclotron frequency. Thus, they can modify the plasma behavior 
in different ways. In particular, charged dust particles introduce a cutoff (resonance) in 
the Alfven wave spectrum at the dust cyclotron frequency. For the ions, this resonance 
occurs in a narrow range of higher frequencies, which are unimportant in the systems under 
consideration. On the other hand, the dust cyclotron resonance occurs at low frequencies 
and it can be an important damping mechanism for the waves. 

If a distribution of grain sizes is considered, then we obtain a band of resonance fre- 
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quencies from the minimum dust-cyclotron frequency (^min to the maximum dust-cyclotron 
frequency cUmax- According to Mathis et al. (1977), we can describe the distribution of grain 
sizes in the interstellar medium by a power law (MRN distribution): 

f(n)dn^Cpn-pdn, (i) 

where Cp — (p — 1)/(1 — a^^), — a/omin is a dimensionless radius, a is the grain ra- 
dius, — Omax/omin, ^ud Omin and Omax ^rc the minimum and the maximum dust radii, 
respectively. Observationally, the parameter p depends on the dust constituent and on the 
environment. In this work we consider p = 4, as used by Cramer et al. (2002). In this case, 
we have C4 = 3/(1 — a~^). 

The wave propagation in dusty plasma is modified and some new and interesting effects 
take place (Wardle & Ng 1999; Cramer 2001). The dispersion relation of the Alfven waves, 
with frequencies smaller than the ion cyclotron frequency, considering dust particles with 
constant charges in a neutral and cold dusty plasma, is given by (Cramer et al. 2002) 

ui±U2, (2) 

where: 

^ u'nj 0^ r™ f{n)dn 

and 

u^^i cu^nao T"^ nf{n)dn 



nf{n)dn 



Here s — C4ln(a^), uj is the angular wave frequency, Qj is the ion cyclotron frequency, 
^do = ?min-B/ (mminc) = (^max IS the maximum dust cyclotron frequency (i.e., for the minimum 
dust grain with mass mmin and charge q'min), and VAa = B / [\/ATipa) is the Alfven speed of 
the a species. Here a = i and d for ions and dust, respectively, B is the magnetic field 
intensity, is the a species density, and c is the speed of light. In order to calculate the 
grain mass, we assume spherical silicates grains with p^n = 3.3 g cm~^. 

Dust grains are usually negatively charged, with many electrons collected on each par- 
ticulate. The mean grain particle charge can be obtained for each dust radius by considering 
the charge current equilibrium over the dust surface. The equilibrium equation, if we consider 
a Maxwellian distribution of velocities, is given by (Vladimirov 1997) 

LO^^ V aksT J vt, [ aksT J 
The plasma frequency, a;^^, for electrons {j3 = e) or ions (/3 = i) is given by: 

lATippql 

^v. = ^^ (6) 
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and the thermal velocity, vxg, is given by: 



'"Tg = \ , (7 



where is the Boltzmann constant. 



The left-hand polarized wave (minus sign in eq. [2] ) interacts with ions and is not affected 
by the negatively charged dust particles resonance. The right-hand polarized wave (plus sign 
in eq. [2]) is the mode damped by the dust resonance. In the case of (^min < uj < cuuiax, the 
integral over the particle's radius has singularities, whose residues give the complex part of 
the wave number (/c^ = kn + ikc) and that leads to the dust-cyclotron damping of the waves, 
with damping length: 

kcioj) 

We adopt a particle radius distribution from = 5.0 x 10~^ cm to flmax = 2.5 x 10^^ cm 
(Rodgers & Glassgold 1991). In order to calculate pd, we adopt a gas-to-dust ratio given by: 

Rgd^-^. (9) 
Pd 

where p is the gas density. 

In this work, we do not develop a self-consistent model for the grain nucleation. Accord- 
ing to the model of Dominik et al. (1989) for dust-driven winds around carbonated stars, 
near the critical point (~ 1.5 tq) there is an intense grain formation. For r < 1.5 Tq, grains 
are only formed locally, and beyond 2.0 Tq, dust is already condensed into an envelope around 
the star. On the basis of this scenario, wc adopt the following situation. For r > 2.0 tq, 
Rgd ~ 200 (Knapp 1985; Heras & Hony 2005), which implies that the dust envelope has 
already been formed. Near the wind base, Rgd 3> 200, which imphes the existence of a 
negligible dust quantity. Between the wind base and ~ 2.0 Tq, Rgd decreases exponentially, 
which indicates that the grains are being formed until Rgd reaches its minimum value of 
~ 200. 



2.2. Magnetic Field Geometry 

It is beyond the scope of this paper to discuss magnetic field generation in late-type 
stars. However, many works studying field generation have been presented in the literature 
(Blackman et al. 2001; Soker & Zoabi 2002; Dorch 2004). These works are motivated by 
observations of, for instance, maser polarization (Vlemmings et al. 2002) or X-ray emission 
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(Ayres et al. 2003) from late-type stars, which are taken as an indicative of the presence of 
magnetic fields around them. Hence, jointly with perturbations (e.g., convective motions), 
the basic ingredient for Alfven wave generation is present in the atmospheres of late-type 
stars. 

On the basis of this assumption, Jatenco-Pereira & Opher (1989) suggested a model for 
mass loss in late-type stars driven by an outward flux of damped Alfven waves that took into 
account our knowledge of coronal holes in the Sun. It is known that the coronal holes are the 
source of the high-speed solar wind streams at the Earth's orbit (e.g., Peter & Judge 1999). 
The role of Alfven waves in accelerating the fast solar wind and the heating of plasma in the 
open magnetic structures of the solar corona has been discussed since the pioneering work 
of Parker (1965) (e.g., Rosenthal et. al 2002; Suzuki & Inutsuka 2005, among others). The 
area expansion of a solar coronal hole is 5 — 13 times greater than that for a radial expansion 
(Lie-Svendsen et al. 2002; Esser et al. 2005; Tu et al. 2005). For simplicity, we assume that 
in the wind of a typical K5 supergiant star, the non-radial expansion factor is 

where flo and fl are the solid angles at the stellar surface and at the edge of the coronal hole, 

respectively. Hence, in order to reproduce the coronal holes, we use a diverging geometry 
suggested by Kuin & Hcarn (1982) for the magnetic field. The cross section of a flux tube 
at a radial distance r is given by 



Air) 



A{ro){r /tq) if r < 



Mro)irt/ro)'^{r/rtf if r > , 

where is a parameter that determines the divergence of the geometry up to the transition 
radius r^, as shown in Figure 1. For a given S, rt is easily obtained: 

Mrt)M fny-' 

Oo A(ro)/r2 \roJ ' ^ ^ 

Hence, we have 

r, = FV(^-2)ro. (12) 

Then, when we consider this geometry, conservation of magnetic flux yields a magnetic field 
intensity B{r) oc A{r)~^. 



2.3. The Wind Equations 

The mass, momentum, and energy equations are the equations responsible in driving 
the wind: 
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The mass equation It expresses the continuity of matter. For a flow velocity u, we have: 

puA{r) = C . (13) 

The constant C in equation (13) is computed at the steUar surface: C — poUoA{ro), where 
the subscript "0" indicates that the parameter is being evaluated at r = Tq. 



The momentum equation In this model, we assume that the forces acting on the wind 
particles are the gravitational force, the gas pressure gradient, and the wave magnetic pres- 
sure gradient. Hence, if we assume a steady flow, the equation of motion can then be written 
as: 

^du _ ^GM^ dP 1 de 
dr dr 2 ' 

where P — pkBT/m is the gas pressure, m is the average mass per particle, e is the energy 
density of the Alfven waves, and GM^jr^ is the gravitational acceleration. The last term 
in equation (14) is the wave magnetic pressure gradient. As dP/dr and de/dr are both 
negative quantities, we can see that the gas and the wave pressure gradients act in the 
opposite direction to the gravitational attraction. 



The energy equation In order to compute the temperature profile of the wind, one needs 
to evaluate the energy equation. In the present model, the wind temperature is determined 
through the balance between wave heating, adiabatic expansion, and radiative cooling. If we 
neglect conduction and assume a perfect gas, we write the energy equation as (Hartmann et 
al. 1982) 

where is the escape velocity. The term {u/'i)de/dr is the rate at which the waves do work 
on the gas. The quantity Q is the wave heating rate, i.e., the rate at which the gas is being 
heated due to the dissipation of wave energy, and Pr is the radiative cooling rate, both in 
units of ergs cm~^ s~^. 

The radiative cooling rate is given by: 

Pr^ AneUH , (16) 

where rie is the electron density, uh is the hydrogen density, and A is the radiative loss 
function. Here, we adopt the value of A given by Schmutzler & Tscharnuter (1993) and 
calculate Ue as did Hartmann & MacGregor (1980). 
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In order to evaluate the wave heating rate Q, one needs to know how the wave energy 
is being dissipated. In our model, we assume that Alfven waves are damped by the dust- 
cyclotron damping mechanism, which has a damping length L{uj) given by equation (8). It 
is importante to note that each wave frequency is damped by a different amount (sec §3 for 
a discussion). Hence, the wave heating rate must be different for each wave frequency. The 
wave heating rate that enters equation (15) is, then, calculated as follows: 

g = / Q{uj)duj, (17) 

where Q{uj) is the wave heating rate per frequency and it is directly related to the wave 
energy density per unit frequency e{uj) and the damping length L{uj) associated with it. 
Hence, we have (HoUweg 1973) 

Q{u^)^'^^{u + va). (18) 

Similarly, the wave energy density e needs to be calculated over the frequency band. 
Hence, we have: 

e= / e{uj)duj (19) 



and, for each frequency, it is calculated as (Jatcnco-Pereira & Opher 1989): 

Mo /i + MoV r 1 , , 

t{u) = eo{uj)— exp / --rr^dr , (20) 



M \1 + M J L{u;) 

where M — u/va is the Alfven Mach number with va — (B/y/ATrp). 

Tu et al. (1989) observed an Alfven wave spectrum propagating in the solar wind and 
inferred for the low-frequency range that e(a;) oc u!~^, with /3 ~ 0.6. Here, we assume that 
the initial energy density of Alfven waves propagating in the winds of cool supergiant stars 
behaves similarly: 

eo(^) = eoiuo) ( —] , (21) 



where eo(cc'o) is the initial wave energy density at any frequency ujq. The initial energy density 
and the initial wave flux evaluated at cu are related by (Jatenco-Pereira & Opher 1989): 

(PAoiuj) = eo{uj)vAo (l + ^^o) • (22) 
In our model, we adopt an initial wave flux as estimated by Hartmann & MacGregor (1980): 

0Ao = / 0Ao(^) duj ~ 10^ erg cm'^ s"^ . (23) 
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Prom equations (14), (15), and (20) we can write the temperature variation as: 

r{Q - Pr) 



dT _2T 
dr 3 r 

and the velocity variation as: 

1 du 

u dr 
_ Z 
r 

3 m 



pu{kBT /m) 



-{z+^y^ 

u dr 



u 



3 m 

3 



Zpu 



+ 



2Zp 



e /1 + 3M 
Ip \1 + M 

e /1 + 3M 
4^ V 1 + M 

L[u!) Z 



+ 



+ 



where we define Z as: 



S 
2 



(24) 



(25) 



if r < rt 
if r > rt . 

Hence, in order to obtain the velocity and the temperature profiles of the wind, we then 
solve equations (13), (24), and (25), simultaneously, for the wind of a typical K5 supergiant 
star. 



2.4. Radiation Pressure 

As in Paper I, we include the radiation pressure into the wind equation using an effective 
escape velocity: 

i^?^^(i-r). (26) 

where F is the ratio between the force per unit volume exerted on rid grains of radius a situ- 
ated at a distance r from a star of luminosity [ridTra^QdL^/ {4:Trcr^)] and the gravitational 
force per unit volume {pCM^./r'^): 

^ UdTra^OdL^ 

p AttcGM^ ' ^ ' 

Here Qd is the mean extinction efficiency factor of the grain. In order to analyze how effective 
the radiation pressure is on the grains, we make a simple estimate. For silicate grains with 
a ~ 5 X 10~^ cm in the atmosphere of a star with an effective temperature of 3500 K, we 
have Qd ~ 10"^ (Oilman 1974). We can write: rid/ p — pd/ip^d) — (J^gdT^d)'^- Then, if 
we assume that TZgd ^ 200, we obtain F ~ 0.03; that is, the radiation pressure modifies the 
square of the escape velocity by only ~ 3%. This is not a significant modification. 
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in Paper I, larger grains were considered. As Qd is more important for large grains, it 
was shown in Paper I that radiation pressure leads to significant higher terminal velocities, 
although it can not initiate the wind. It is observed that the wind forms before the grain 
formation point (Carpenter et al. 1995), which indicates that another mechanism is taking 
place in order to drive the wind. 

In the present model, we observe only a negligible difference in the terminal velocity if 
we consider the radiation pressure or not. This happens because we are adopting the radius 
distribution given by equation (1), which causes our model to have many more small grains 
than large ones. 

3. RESULTS AND DISCUSSION 

For the solution of the set of equations (13), (24), and (25), we adopted the set of initial 
conditions presented in Table 1 for the wind of a typical K5 supergiant star. We solve these 
equations until r = 300 tq, at which point the wind has already reached its terminal velocity. 
Some parameters adopted here, such as the stellar mass, the stellar radius, the initial density, 
and the initial magnetic field, all have of the same values as the ones adopted by Hartmann 
& MacGregor (1980) for their model 6. The parameters calculated for the dust grains are 
shown in Table 2. 

The momentum equation (eq. [25]) of the wind has a critical point when both the 
numerator and denominator go to zero. The requirement that the wind velocity should 
increase through the critical point determines the initial velocity at the base of the wind, 
i.e., du{r — Tq) /dr > 0. Hence, to find the critical solution for a given set of initial conditions, 
we iterate the initial velocity uq until the solution passes through the critical point. Together 
with the initial density, the initial velocity determines the mass-loss rate. 

To simplify the problem, we calculate cu^in and cu^ax at r = Tq and assume that these 
parameters are constant throughout the wind. The grains charges are evaluated according 
to equation (5). Our results are shown in Figures 2 and 3, where we plot the velocity and the 
temperature profiles that we obtained. In order to have a better view of what is happening 
near the stellar surface, we plot the profiles until r = 30 Tq. 

The results of the present model compares favorably with observations. We expect 
the terminal velocity to be lower than the surface escape velocity. From observations, we 
expect Uoo ~ Veo/2. Our model provides Uoo — 57 km s~^ and M ~ 2.5 x 10~^ Mq yr~^. In 
Figure 2, we can see that the wind is accelerated near the surface, even before the grain 
formation point (~ 2 Tq). 
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As we know from solar corona observations, we expect a rise in the wind temperature 
near the stellar surface. From Figure 3, it can be seen that the rise in the temperature occurs 
very abruptly near the wind base. This heating is explained by the mechanical dissipation of 
Alfven waves. At larger distances (r > 30 Tq), the flow expands adiabatically, hence, in the 
absence of strong heating, the temperature tends to fall monotonically with the adiabatic 
exponent 4/3. 

In order to analyze the wave dissipation, we plot in Figure 4 the energy density of 
the waves at different distances from the star. At r = Tq, we have a power-law spectrum 
(eq. [21]). As wc get farther from the wind base, the energy density decreases due to the ffow 
expansion and wave damping. In Figure 4, we note that this decline is larger for larger wave 
frequencies, indicating a more efficient damping for waves with higher frequencies. This is 
due to the fact that we have a larger number of small-radius grains (see eq. [1]). As the 
smaller grains have larger cyclotron frequencies, they interact with the higher frequency 
waves, causing a larger damping at this spectral range. 

Jatenco-Pereira & Opher (1989) studied the mass loss in late-type supergiants using 
three different damping mechanisms: nonlinear damping, resonant surface damping, and 
turbulent damping. They showed that these damping mechanisms have similar behaviors 
in the wind. In Figure 5, we compare the wave flux (J)a calculated by Jatenco-Pereira & 
Opher (1989), damped by the resonant surface damping, and the frequency- integrated flux 
calculated by the present model; that is, damped by the dust-cyclotron mechanism. It can 
be seen that the damping of the waves by the dust-cyclotron mechanism is not as efficient as 
any of the others mechanisms studied by Jatenco-Pereira & Opher (1989): at ~ 4 ro the wave 
flux has droped by more than a factor of 10^ in the model of Jatenco-Pereira & Opher (1989). 
Also, the dust-cyclotron damping mechanism is not as strong as the exponential decay that 
was simulated in Paper I. Nevertheless, our model can reproduce observational data such as 
the mass loss rate and the terminal velocity for the wind of a typical K5 supergiant star. 



4. CONCLUSIONS 

In this work, we consider an outward-directed flux of damped Alfven waves as an accel- 
eration mechanism for driving late- type stellar winds. As these winds present great amounts 
of dust particles and as the dust can modify the wave modes, we study here how the dust can 

affect the propagation of Alfven waves, taking into account a speciflc damping mechanism, 
dust-cyclotron damping, which can be understood as follows. A charged grain immersed in 
a magnetized plasma possesses a cyclotron frequency associated with it. If Alfven waves are 
present in this plasma, their propagation near the dust-cyclotron frequency can be affected. 
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As different grain sizes carry different charges, there will be a cyclotron frequency associ- 
ated with each grain size. In this sense, a dust size distribution has several dust-cyclotron 
frequencies, which will affect a broad band of Alfven wave frequencies. 

We also calculated the radiation pressure exerted on grains. We included this force as 
another acceleration mechanism of the wind. As we adopted great amouns of small grains, 
we observed only a negligible difference in the wind velocity profile from what we would have 
obtained if we did not consider radiation pressure. 

Our model also assumes a diverging geometry for the magnetic field lines. This assump- 
tion is based on the observations of coronal holes in the Sun, which present a superradial 
expansion of magnetic field lines. Thus, the mass, momentum, and energy equations are 
obtained and then solved in a self-consistent approach. 

Our results of wind velocity and temperature profiles for a typical K5 supergiant star 
shows compatibility with observations. We also show that, considering the presence of 
charged dust particles, the wave flux is less damped due to dust-cyclotron damping than 
it would be if we considered some other damping mechanisms studied in the literature, such 
as nonlinear damping, resonant surface damping, and turbulent damping. 

The authors would like to thank the Brazilian agencies FAPESP (under grant 02/10846-0) 
and CNPq (under grant 304523/90-9) for financial support. 
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Fig. 1. — Geometry used in order to reproduce coronal holes (not to scale). 
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7000 




Fig. 3. — Wind temperature profile obtained for a typical K5 supergiant star. 
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Fig. 4. — Energy density of Alfvcn waves at different positions: at 1 tq (solid curve), at 
1.34 ro (dotted curve), at 2 tq (dashed curve), at 2.34 tq (dash-dotted curve), and at 5 Tq 
(dash-double-dotted curve). 




Fig. 5. — Comparison between Alfven waves fluxes for tliis model (solid curve) and for the 
model of Jatenco-Pereira & Opher (1989) (dash-dotted curve). 
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Table 1: Initial parameters adopted for the wind of a typical K5 supergiant star. 



Parameter 


Value 


Unity 


To 


400 






16 


Mo 




5.5 X 10^ 


erg cm~^ s"^ 


Po 


1.07 X 10-13 


g cm~3 




1.5 X 10^ 


cm s"-*^ 


Bo 


10 


G 


To 


3500 


K 


S 


4.0 




n 


3.16 


To 



Table 2: Maximum and minimum calculated initial dust charges and cyclotron frequencies 
related to the minimum grain radius (omin = 5.0 x 10"^ cm) and the maximum grains radius 
(omax = 2.5 X 10-^ cm). 



Parameter 


Value 


Unity 


?min 


2.6 


e 




131.2 


e 


^min 


9.7 X 10"^ 


rad s"! 


^max 


2.4 X 10-1 


rad s-i 



